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NOMENCLATURE 


A  vector  is  denoted  with  an  over-arrow  (e.g.,  the  position  vector  the  column  matrix  rep¬ 

resentation  of  a  vector,  or  any  other  type  of  column  matrix,  is  represented  with  an  underbar  (e,g.,  v). 
Points  in  the  system  are  designated  by  lowercase  letters,  so  the  position  vector  in  the  above  example 
is  read  “the  position  from  point  a  to  point  b.”  Coordinate  frames  are  designated  by  capital  letters  (e.g., 
XX  is  the  X  unit  vector  in  the  A  system).  Coordinate-frame  rotation  rates  are  represented  as  uJjr//, 
which  is  read  as  “the  rotation  rate  of  frame  F  relative  to  frame  I.”  Rotation  vectors  are  always  resolved 
in  the  frame  of  the  first  subscript.  Nondimensional  quantities  are  represented  with  an  overbar  and  are 
nondimensionalized  by  the  blade  mass,  rotor  radius,  and  nominal  rotor  rotational  rate  .The  sine  and 
cosine  functions  are  abbreviated  to  Sq  and  Cq,  where  the  subscript  is  the  argument.  In  the  following 
list,  the  abbreviation  DOF  is  used  for  degrees  of  freedom. 
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SUMMARY 


A  mathematical  model  of  a  helicopter  system  with  a  single  main  rotor  that  includes  rigid,  hinge- 
restrained  rotor  blades  with  flap,  lag,  and  torsion  degrees  of  freedom  is  described.  The  model  allows 
several  hinge  sequences  and  two  offsets  in  the  hinges.  Quasi-steady  Greenberg  theoi7  is  used  to 
calculate  the  blade-section  aerodynamic  forces,  and  inflow  effects  are  accounted  for  by  using  a  three- 
state  nonlinear  dynamic  inflow  model.  The  motion  of  the  rigid  fuselage  is  defined  by  six  degrees  of 
freedom,  and  an  optional  rotor  rpm  degree  of  freedom  is  available.  Empennage  surfaces  and  the  tail 
rotor  are  modeled,  and  the  effect  of  main-rotor  downwash  on  these  elements  is  included.  Model  trim, 
linearization,  and  time-integration  operations  are  described  and  can  be  applied  to  a  subset  of  the  model 
in  the  rotating  or  nonrotating  coordinate  frame,  A  preliminary  validation  of  the  model  is  made  by 
comparing  its  results  with  those  of  other  analytical  and  experimental  studies.  This  publication  presents 
the  results  of  research  completed  in  November  1989. 

INTRODUCTION 

Future  requirements  for  helicopter  flight  performance  place  greater  emphasis  on  high-bandwidth 
flight-control  systems  than  current  requirements.  Low-order  rigid-body  models  of  helicopter  systems  can 
be  used  successfully  to  design  feedback  control  systems,  but  these  models  impose  bandwidth  limitations 
on  the  controller,  making  them  unsuitable  for  the  design  and  evaluation  of  high-bandwidth  systems 
(ref.  1).  To  meet  the  new  requirements,  higher  order  helicopter  models  must  be  used  in  the  design  of 
the  flight-control  systems.  The  higher  order  effects  of  primary  interest  are  those  of  the  rotor  dynamics 
and  the  rotor  inflow.  In  addition  to  capturing  the  blade-flap  and  lag  dynamics,  the  effect  of  torsional 
dynamics  on  the  control  design  problem  needs  to  be  understood.  This  last  item  is  relevent  to  rotor 
designs  using  hingeless  and  bearingless  rotors  and  to  advanced  concepts,  such  as  servo-flap  control, 
which  require  torsionally  soft  rotor  blades  (ref.  2). 

This  report  documents  a  higher  order  helicopter  mathematical  model  created  to  fill  these  needs.  The 
helicopter  model  is  described  in  detail  in  the  next  section.  The  rotor  consists  of  rigid  blades  with  flap, 
lag,  torsion,  and  pitch  motions;  and  with  linear  hinge  springs  and  dampers.  This  gives  the  model  the 
rotor  dynamics  of  interest  and  allows  for  approximate  modeling  of  hingeless  rotor  systems.  The  flap  (f), 
lag  (1),  torsion  (t),  and  pitch  (p)  hinge  motions  can  have  the  order  Apt,  fplt,  pflt,  Ifpt,  Ipft,  or  plft  (only 
Apt  and  Ifpt  are  currently  available),  and  two  hinge  offsets  are  available  to  model  various  articulated 
configurations.  The  blade-section  aerodynamics  are  linear,  and  the  quasi-steady  Greenberg  model  is  used 
(ref.  3).  Unsteady  inflow  effects  are  included  using  the  three-state  nonlinear  Pitt/Peters  dynamic  inflow 
model  (ref.  4).  An  rpm  degree  of  freedom  (DOF)  of  the  main  rotor  is  also  available,  making  analysis 
of  engine/rotor  interaction  possible.  The  rigid  fuselage  model  has  six  degrees  of  freedom  and  includes 
an  equivalent-drag-area  aerodynamic  model,  a  tail-rotor  model,  and  an  empennage-surface  model.  The 
empennage-surface  models  have  main-rotor  downwash  effects  that  are  calculated  from  the  flat-wake 
model  described  in  chapter  2  of  reference  5,  The  tail  rotor  is  the  quasi-static  flapping  model  described 
in  reference  6.  For  greater  flexibility,  the  spring/damper  model,  the  blade-section  aerodynamic  model, 
and  the  fu.selage  aerodynamic  model  can  be  replaced  with  user-defined  models. 

Following  the  description  of  the  model  is  the  Operations  section,  which  details  the  operations  that 
are  performed  on  the  model.  There  are  three  operations:  trim,  linearization,  and  time-integration.  Trim 


places  the  helicopter  system  in  some  user-defined  flight  condition  by  satisfying  the  system  equations  of 
motion.  After  trimming,  the  model  may  be  linearized  about  the  trim  condition,  or  the  model  may  be 
initialized  in  the  trim  condition  and  then  time-integrated.  Linearization  uses  nonrotating  rotor  coordinates 
through  a  multiblade  coordinate  transformation  (MBCT)  (ref.  7).  In  forward  flight,  the  periodic  system 
matrices  can  be  extracted  from  this  linearization  process.  The  time-integration  uses  rotating  blade 
coordinates  and  allows  for  arbitrary  inputs  to  the  model  controls.  Additionally,  the  integration  process 
has  been  cast  in  a  general  block  form,  allowing  the  introduction  of  user-specified  actuator,  sensor, 
feedback,  and  precompensator  dynamics  to  the  system. 

In  the  Results  section,  a  preliminary  validation  of  the  helicopter  model  is  made  by  comparing 
results  with  those  from  other  analytical  and  experimental  studies.  The  rotor  model  results  are  compared 
with  trimmed  rotor-response  results,  from  reference  8,  of  an  isolated  rotor.  Rotor/fuselage  results  are 
compared  with  the  experimental  frequency  and  damping  data  presented  in  reference  9.  The  full  helicopter 
model  is  validated  by  comparing  trim  and  time-history  data  extracted  from  flight  data  on  a  UH-60A 
helicopter  (ref.  10).  Additionally,  a  comparison  with  the  quasi-static  eight-state  linearized  model  from 
the  GEN  HEL  program  (refs.  11,12)  is  made. 

HELICOPTER  EQUATIONS  OF  MOTION 
Fuselage 

The  fuselage  is  modeled  as  a  rigid  body  whose  motion  is  defined  by  the  six  degrees  of  freedom 
shown  in  figure  1.  The  fuselage  coordinate  frame  F  is  body  fixed,  with  its  origin  at  point  m.  The 
location  of  point  m  is  chosen  by  the  user.  The  body  translational  motion  is  defined  by  the  components 
of  the  inertial  position  vector  rj,n*  and  body  rotation  is  defined  by  the  Euler  angles.  Other  relevant  points 
on  the  body  are  shown  in  figure  2.  The  positions  of  these  points  are  chosen  by  the  user  in  the  body-fixed 
coordinate  frame  defined  by  the  station,  the  buttline,  and  the  waterline.  The  forces  and  moments  at 
m  result  from  mounting  restraints,  which  can  be  set  to  simulate  any  mounting  condition,  such  as  in  a 
wind-tunnel  test.  Point  c  is  the  center  of  mass  of  the  fuselage,  where  inertial  and  gravitational  forces 
are  resolved,  and  point  a  is  where  fuselage  aerodynamic  forces  are  applied.  The  n.,  empennage  surfaces 
produce  resultant  forces  at  the  sj  points,  and  the  resultant  tail-rotor  force  is  at  point  r.  There  are 
points  designated  as  yn,  where  accelerations  are  output  measurements.  Point  h  is  the  attachment  point 
of  the  main-rotor  hub,  which  can  be  set  at  tilt  angles  about  point  t.  These  tilt  angles  are  defined  by 
the  transformation,  shown  in  figure  3,  which  tips  the  hub  coordinate  frame  H  by  constant  pitch  and  roll 
angles,  9s  and  0^. 

The  equations  of  motion  for  the  fuselage  are  formed  by  summing  the  forces  created  by  the  main 
rotor,  fuselage  inertia  and  gravity,  fuselage  aerodynamics,  the  fuselage  mount,  the  tail  rotor,  and  the  tail 
surfaces.  This  gives  the  vector  equation 

Qh  +  Qig  +  Qa  +  Qm  +  Qr  +  ^  Qsj  =  o  ( i ) 

j 

Summing  the  respective  moments  about  point  rn  gives 

+  ^IG  +■  -^.4  "i'  +  -^T?  +  -^5  “0  (2) 
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where 
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The  final  fuselage  equations  are  resolved  such  that  they  correspond  to  the  fuselage  generalized  coordi¬ 
nates.  That  is,  if  a  Lagrangian  approach  had  been  used,  the  same  fuselage  equations  would  have  resulted. 
This  resolution  of  the  equations  facilitates  the  extraction  of  the  system  mass  matrix,  which  is  discussed 
in  a  subsequent  section.  The  definition  of  the  equations  and  their  location  in  the  fuselage-equation 
vector  are 


moments  about  Stf 
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(9) 

feqi'i)  =  5Znion^ents  about 

(10) 

/e?(4)  =  direction 
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/eg(5)  =  forces  in  y/ direction 

(12) 

/eg  (6)  =  5Z  forces  in  1/ direction 

(13) 

Using  the  position  vectors  and  coordinate  frames  as  defined,  the  velocities  and  accelerations  can 
be  determined.  At  the  main  rotor  hub,  the  angular  velocity  and  acceleration  are  given  by 


<^H/I  =  ‘^H/F  +  '^F/I 

(14) 
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(15) 
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Taking  advantage  of  the  constant  components  in  the  position  vectors,  the  velocities  at  selected  points 
are 

ru  =  Tim  +  X  Tmt  (17) 
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It  It 
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(27) 

The  forces  and  moments  of  the  fuselage  are  discussed  separately,  beginning  with  the  inertial  and 
gravity  effects.  At  the  fuselage  center  of  gravity,  the  inertial  force  and  moment  vectors  are 

QiG  -  (“He  +  qZi)  (28) 

.7/  =  —  Ip  -  Fjl  —  I^F/I  ^  i^F  '  '^F/l) 
where  the  inertia  tensor  is  defined  in  the  F  system. 


The  fuselage  aerodynamic  effects  are  based  in  the  wind-axis  coordinate  frame  W  defined  in  figure  4. 
The  velocity  of  the  air  relative  to  point  a  is 

Va  =  -  ha  (30) 

which  also  defines  the  direction  of  the  xyp  component  of  the  wind-axis  frame.  The  angles  c\f  and  f3f 
are  defined  from  the  relative  velocity  as 


af  =  tan 
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=  tan  ' 
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The  fuselage  aerodynamic  force  and  moment  are 

Qa  =  +  yyw  +  =  ^pWdW'^Sf  (Cdj£\v  +  (^VjfJw  +  (33) 

Ja  -  -Ri'w  +  Mijw  -  Nzw  =  {-CrjXw  +  -  C^jZ\y)  (34) 
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which  can  be  converted  to  the  F  frame  via  the  transformation  Tp\y,  which  is  defined  in  figure  4.  The 
force  and  moment  coefficients  are  functions  of  of  and  j3f,  but  the  default  fuselage-aerodynamics  sets 
all  force  coefficients  to  zero  except  which  is  set  to  unity,  making  the  fuselage  reference  area  Sf 
the  equivalent  drag  area. 

The  fuselage  mount  moments  are  modeled  as  linear  nonorthogonal  springs,  and  the  mount  forces 
are  treated  as  given  quantities.  These  forces  can  be  used  during  trim  to  satisfy  a  constraint  condition, 
such  as  setting  a  certain  amount  of  thrust.  The  parameters  defining  the  mount  effects  can  be  set  to  zero 
for  free-flight  analysis.  The  force  and  moment  expressions  for  the  mount  are  given  by 

Lm  =  {Mx  -  MzSg^')  Xp  -j-  [MijCg^  +  MzSg^COy)  Vp  +  (^-AfySg,^  +  zp  (35) 

where 

Mx  =  -Kp.^  [6x  -  Ox^ )  -  Cp.^^  ix  (36) 

A4  =  -Kf^  («„  -  e„„)  -  Cf/„  (37) 

M,  =  -Kf,  (fl,  -  «,„)  -  (3R) 

Qm  -  +  FyOj  -1-  F^zj  (39) 

The  measurement  outputs  from  the  fuselage  are  acceleration,  velocity,  and  angular  rates.  The 
accelerometer  outputs  are  expressed  as 


;v(i)  =  -  gzij- 


j  =  1,2,.. .  ,7? 


fj  =  fxjXp  +  JyjVp  4-  fzjh^ 


(40) 

(41) 


where  the  direction  of  the  measurement  is  selected  through  the  components  of  fj.  The  velocity  of  the 
fuselage  center  of  mass  is  resolved  in  the  fuselage  fixed  frame  and  is  available  as  a  measurement,  as 
are  the  fuselage  angular  rates.  The  components  of  the  velocity  and  angular  rate  vectors  are  stored  in  y 
according  to 


11 

+ 

(42) 

y  +  2)  =  He  •  VF 

(43) 

y  +  3)  =  fic  ■  zp 

(44) 

11 

+ 

(45) 

11 

''iFT' 

+ 

(46) 

11 

(47) 
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Main  Rotor 


Fixed  to  the  main  rotor  are  coordinate  frames,  one  for  each  blade.  Figure  5  shows  the  rotor 
coordinate  frame  R  for  the  kth  blade,  in  which  several  additional  blade  coordinate  frames  are  defined. 
The  additional  frames  are  shown  in  figure  6.  The  first  hinge  is  offset  a  distance  e  from  the  hub  center, 
and  point  c  corresponds  to  the  origin  of  the  link  coordinate  frame  L.  The  second  hinge  is  offset  a 
distance  /  from  point  f  .  The  origin  of  blade  coordinate  frame  B  is  at  point  /  and  is  fixed  to  the  blade. 
The  significant  points  on  the  blade  are  b,  at  the  blade  center  of  gravity,  and  .r,  at  a  distant  .v  along 
the  d'n  axis.  Also  on  the  blade  are  points  designated  tcn.  which  locate  the  blade  at  erometer 
measurements. 

The  hinge  sequences  that  define  the  orientation  of  the  blade  coordinate  frames  can  be  divided  into 
two  groups  .  The  first  group  is  flap-lag-pitch-torsion  (Apt),  flap-pitch-lag-torsion  (fplt),  and  pitch-flap- 
lag-torsion  (pflt),  in  all  of  which  the  flap  motion  precedes  the  lag.  The  second  group  places  lag  before 
flap  in  lag-flap-pitch-torsion  (Ifpt),  lag-pitch-flap-torsion  (Ipft),  and  pitch-lag-flap-torsion  (plft)  hinge 
sequences.  Only  the  Apt  and  Ifpt  sequences  are  available;  they  are  defined  in  figures  7  and  8.  An 
intermediate  coordinate  frame  P,  which  is  not  shown  in  figure  6,  is  included  in  these  figures.  Frame 
P  allows  the  introduction  of  hinge  effects  such  as  precone  and  (^3.  The  angle  (Sp  in  figures  7  and  8  is 
the  precone  angle,  and  Cp  is  a  cant  in  the  flap  hinge  that  geometrically  introduces  the  ^'3  effect.  By 
definition,  the  variable  6^  is  not  an  Euler  angle  (ref.  13),  but  it  is  related  to  the  cant  Euler  angle  by 
tan  6^  =  -  tan  (p  cos  j3p  for  the  Apt  system  and  tan  <53  =  -  tan  +  Cfc)  cos  f3p  for  the  Ifpt  system.  The 

remaining  hinge  sequences  are  defined  such  that  for  the  flap-before-lag  sequences, rotates  about  the 
coordinates  axis  -j/p;  Cfct  about  and  about  xg.  For  the  lag-before-flap  sequences,  rotates 
about  the  axis  £p;  about  -{/£,;  and  about  xg-  The  sequence  of  rotations  of  the  blade  angles  is 
summarized  in  table  1.  The  items  in  the  Apt  row  and  R  to  P  column  are  read  as  follows:  starting  in 
frame  R.  rotate  through  an  angle  of  /ip  about  the  negative  yg  axis;  then  rotate  through  an  angle  of  Cp 
about  the  axis,  to  the  P  system.  Each  element  of  the  table  is  read  similarly. 

The  forces  on  the  main  rotor  are  shown  in  figure  9  with  only  the  A'th  blade  depicted.  A  rigid  body 
representing  the  rotor  hub  mass  contributes  only  inertia  to  the  system.  The  link  between  points  e  and  / 
is  massless,  contributing  no  inertial  forces  to  the  system.  Each  blade  contributes  inertial,  gravitational. 


and  aerodynamic  forces  and  moments,  which  are  related  to  the  loads  at  hinges  /  and  e  by 

Ff  =  -  {Fjg  +  Fa)  =  -Ft  (48) 

Mf  =  -  (Mig  +  Ma)  =  -Mt  (49) 

Fe  =  F/  (50) 

M.=:i\if  +  r,fxFf  (51) 

The  vectors  Fg  and  A/j  are  the  sums  of  the  blade  inertial,  gravitational,  and  aerodynamic  forces  and 
moments  at  point  /.  The  hub  forces  and  moments  of  equations  (1)  and  (2)  for  the  main  rotor  are  then 

Qh  =  F/c  +  i:  (Fr)  (52) 
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Jh  -  [(f/ie  +  fg/)  X  fV  +  -^r] , 

(53) 

k=l 

This  gives  the  rotor  equation 

feqil)  =  Tq  +  Jii  •  ZR  =  moment  about  f/j 

(54) 

where  Tq  is  the  main-rotor  torque.  For  the  Apt,  fplt,  and  pflt  systems,  the  flap,  lag,  and  torsion  equations 

are 

Ttq{h,  1)  =  -  [Me  -f  My  +  fef  Ft]'  VP 

(55) 

req{k,2)  =  [m_^  +  Mj-]  •  zi 

(56) 

reg(fc,3)=  [M_^-f-M7’]  •  xq 

(57) 

For  the  Ifpt,  Ipft,  and  plft  systems,  these  equations  are 

req{k,  l)  =  -[Mf  +  Kif]  •  VL 

(58) 

rtq{k,  2)  =  [Me  +  Mr  +  X  Ft]  •  zp 

(59) 

r,q{k,3)=[Mf  +  MT]'  ip 

(60) 

These  rotor  blade  equations  have  been  written  such  that  they  are  equal  to  the  fik,  4,  and  (l>k  equations 

that  would  result  from  using  a  Lagrangian  approach.  The  total  aerodynamic  forces  and  moments  at  the 

hub  center  are 

Sf,  =  E  {^a), 

(61) 

k=[ 

+  {ne+ref)  X  Fa]i^ 
k==\ 

(62) 

The  angular  velocities  and  accelerations  are 

’^R/[  = 

(63) 

upil  =  iX'p/fi  +  Up/i 

(64) 

-^'LIP  +‘^P// 

(65) 

‘-I 

-h 

II 

(66) 

^Rll  =  Rjl  -  <^'rIH  +  ^'Rll  ^  ^RIH  +  ^'HII 

(67) 

IXpjl  =  =  ix'p/p  +  ^'p/I  X  tX'p^p  +  LCpjl 

(68) 

^'Ljl  =  Lj I  =  L/P  ^  ^'L/P  +‘^P// 

(69) 

^'B/I  =  =  ‘^’'b/L  •+■  ^Bjl  ^  ‘^B/Z,  +  ^'L/l 

(70) 
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For  the  blade  aerodynamic  calculations,  the  velocity  of  point  x  results  from 

^te  “  ^ih  '^R/I  ^  ^/le  ^ ) 

^if  =  'rie  +  ‘^LII  (72) 

'^'ix  ~  "I"  ^B/I  ^  ^/x  (72) 

For  the  mass  matrix  calculation,  the  velocity  of  the  blade  center  of  mass  is 

Tib^fif  +  us^jxffi,  (74) 

The  accelerations  of  the  blade  center  of  mass,  the  points  Wn,  and  the  point  x  along  xq  result  from 

he  =  m  +  uj'fifj  X  f)ie  +  <^R/I  ^  {^R/1  ^  4e)  (75) 

hf  =  he  +  ^  ^ef  +  ^L/I  X  {^L/I  X  r^f)  (76) 

hb  =  f*!/  +  uj'b/I  X  f/6  +  i^B/I  X  {^B/I  X  ^fb)  (77) 

+  ^'b/I  X  +  <^B/I  X  (wj3/;  X  (78) 

hx  =  hf  +  ^'b/I  X  ^/x  +^BII^  ipBIl  X  r/x)  (79) 


The  rotor  hub  center  of  mass  is  at  point  /i,  and  the  rotor  hub  orientation  is  referenced  to  the  fc  =  1 
rotor  coordinate  system.  The  inertial  and  gravitational  effect  of  the  rotor  hub  is 

Pig  =  ^R  {-Uh  +  9^i) 

and  the  inertial  moment  is 

Ri  =  -  [Tr  •  <^'r/i  +  <^r/i  X  [1r  • 

Likewise,  for  each  blade,  the  inertial  and  gravitational  force  and  moment  are 

Pig  ~  ^^B  {-^ib  +  9k) 

^7G  =  -[Tb‘  '^'b/I  +  ‘^S//  X  (Jb  ■  ^B/l)]  +  Pfb  X  ^/G 

where  the  moment  is  about  point  /. 

Blade  hinge-spring  restraining  moments  are  of  two  types.  The  first  type  is  nonorthogonal,  and  for 
the  Apt,  fplt,  and  pflt  sequences  the  moments  are 

“Me  •  yp  =  -Ky(3k  -  Cy^k 
Mf  •  ijr  =  -Kz(^k  -  Cz^k 
Mf  XB  =  -Kx4>k  -  Cx^k 


(80) 

(81) 

(82) 

(83) 
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(84) 

(85) 

(86) 


For  the  Ifpt,  Ipft,  and  plft  sequences,  the  nonorthogonal  moments  are 


(87) 

(88) 
(89) 


-Mf  .  n  =  -Ky/Sk  -  Cyk 

Me  •  ip  =  -A'Kfc  -  C.a 

Mf  •  xp  =  -Kx<f>k  -  Cx^k 


The  second  type  of  spring  restraint  is  orthogonal  and  is  used  to  approximate  hingeless  rotor-blade 
systems.  These  orthogonal  moments  are  valid  only  for  (1)  no  second  offset  (/  =  0),  (2)  no  cant  angle 
((p  =  0),  and  (3)  with  the  pitch  inboard  or  outboard  of  the  flap  and  lag  degrees  of  freedom.  For  the 
Apt  and  pflt  sequences,  the  moments  are 


Mf  =  M,  =  [-A'.  -  Cx  {-CkSp,  +  hC’c.Cfi,)]  xp 

+  [-A'j/  +  <l>kSi^,^)  -  Cy  (-/3fc  -t-  h^Ck)]  VP 

'b  “A'-  (CitC'^fc  +  (l>kC(^^Sj3i^')  -  Cz  (Ck^0it  + 


(90) 


For  the  Ifpt  and  plft  sequences,  the  moments  are 


Mf  =  Me  = 


-Kx  {l3kS^,  +  -  Cx  (4.%  +  hC0,C^,)]  xp 

+  [-A’y  +  ^kCi3k^<:k)  ~  yp 

+  [-As  (Cfc  +  <f>kS0,)  -  Cz  (4  +  zp 


(91) 


The  total  blade  aerodynamic  effects  are  found  by  integrating  the  sectional  aerodynamic  forces  and 
moments  along  the  span  according  to 


-  /*2  - 
=  h<^x 

Jx\ 

Ma  =  [(a;  xp)  X  fA  +  m.4 


(lx 


(92) 

(93) 


where 


Ia  =  fycVC  +  f:c~C 

friA  =  mx^xc 


(94) 

(95) 


The  integrations  are  carried  out  using  Gaussian  quadrature  numerical  integration.  Blade-section  forces 
are  resolved  in  a  chord  coordinate  frame  C  shown  in  figure  10,  which  is  rotated  from  the  blade  coordinate 
frame  B  by  a  geometric  twist.  The  linear  twist  function  is  defined  by 

^T^x)  —  TujxX  -(-  Tujc  (96) 
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The  components  of  the  sectional  forces  and  moments  are  expressed  in  terms  of  the  angular  velocity  and 
acceleration  of  the  chord  section  and  of  the  components  of  the  relative  velocity  and  acceleration  of  the 
air  mass  at  the  origin  of  the  C  frame.  The  latter  two  are 

Vc  =  Vfe  +  (97) 

Ac  =  fG-hx  (98) 

where  the  vectorA^  is  the  induced  air  velocity  whose  derivative  has  been  ignored  in  the  air-mass  ac¬ 
celeration.  The  induced  air  velocity  results  from  the  dynamic  inflow  model,  which  is  described  in  a 
subsequent  section. 


The  sectional  aerodynamic  model  is  based  on  quasi-steady  Greenberg  theory,  which  is  a  Theodorsen 
theory  modified  to  account  for  lead-lag  motions  (refs.  3,14).  The  lift  and  moment  about  the  torsion  axis 
are  given  by 


where 


L  =  Pj\ah 


M  —  Pj\ab 


X'  Q  •+■  P 2^ Aero 

m 

+  2^  “  ^2^  ~  8  2^ Aero 


Q  =  [h  +  Va)  - 
P  =  (jiA-  Vaj  -  Xaja 


(99) 

(100) 


(101) 

(102) 


The  variable  iAero  is  a  switch  to  remove  the  apparent  mass  effects,  which  can  be  ignored  when  the 
order  of  magnitude  of  the  mass  of  the  airfoil  section  is  much  greater  than  the  order  of  magnitude  of  the 
mass  of  the  air  it  occupies  (ref.  14,  chap.  S).  Drag  and  additional  moments  due  to  camber  are  added  by 

D  =  <iQ  + di  oigff  +  ^2  (103) 

^Icam  =  mQ  +  mi  Ogff  (104) 


Figure  11  shows  the  application  of  the  sectional  aerodynamic  model  in  forward  and  reverse  flow  con¬ 
ditions.  The  variables  in  the  Greenberg  theory  are  interpreted  as  in  reference  15,  where  the  quantity 
h  -t-  Va  is  interpreted  as  the  normal  air  velocity  at  the  torsion  axis.  The  derivative  of  this  quan¬ 
tity  is  then  the  derivative  of  the  sc  component  of  the  air  velocity  V'c,  which  can  be  extracted  from 
=  V'c  +  ^  de¬ 


filade  measurements  are  stored  in  w  in  the  following  manner.  For  the  A;th  blade,  the 
acclerometer  outputs  are  given  by 

tD[A;-l-n6(n  -  1)]  =  -f/f/)  •  n  =  1, 2, . . .  ,  n  4^,^;  (105) 

II  II 

+  ^ynVB  +  (106) 
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where  the  direction  of  the  measurement  is  selected  through  the  vector  6n.  Blade  moments  at  the  hinges 
are  also  in  w.  They  are  given  by  the  following  expressions  when  flap  is  before  lag: 


w  [fc  +  =  (Mj  -1-  fef  X  Fj)  •  yp 

(107) 

w  [*  +  fib  +  l)]  =  Mj-  • 

(108) 

w  +  Tib  {^Abid  +  2)  j  =  •  XB 

(109) 

and  by  the  following  expressions  when  lag  is  before  flap: 

w  [i*  +  ribriA^^^]  =  A/j  •  yi 

(110) 

w  [A:  +  Tib  (nAbij  +  01  “  ^ 

(111) 

w  [A;  4-  nb  (nA^,j  +  2)]  =  Mr  •  xp 

(112) 

Inflow 

The  rotor  inflow  is  modeled  using  the  three-state  nonlinear  Pitt/Peters  model  of  reference  4,  which 
is  applied  in  the  actuator-disk  coordinate  frame  A  shown  in  figure  12.  This  frame  is  tipped  to  the  left 
by  the  sine  of  the  flapping  angle  and  tipped  forward  by  the  cosine  of  the  flapping  angle.  The  velocity 
and  acceleration  of  the  air  relative  to  the  hub  center  are 

II 

1 

(113) 

II 

Ci 

1 

(114) 

The  advance  ratio,  the  free-stream  inflow,  and  the  total  average  inflow  can  be  found  by  using  the  velocity 
vectors:  they  are 

m 

(115) 

^  Vh-  ^A 

na 

^  =  ^fs  +  ^0 

(116) 

(117) 

where  the  total  average  inflow  A  is  the  sum  of  the  free-stream  inflow  plus  the  constant  induced-inflow 
component,  which  is  one  state  of  the  dynamic  inflow  model.  The  inflow  vector  for  the  blade  aerodynamic 
calculation,  in  terms  of  the  Inflow  states,  is 

A,'  =  Air^  =  ^Aq  -f-  Ais^50^,  -f-  Aic^Q,^.^  i-l 

(118) 
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where  Aq,  Ai^,,  and  Ai^  are  the  three  inflow  states  that  define  the  inflow  in  the  A  system.  These  inflow 
states  are  governed  by  the  first-order  inflow  equations 


^«g(l)  1 

f  ^  ^ 

^0 

f  1 

1  f  1 

<  Xeq{2) 

^  =  LiMi  < 

♦ 

^Is 

*  +  1 

-i/  Cl 

.  ®fig(3)  J 

* 

^\c  J 

i  J 

1  1  Cm  J 

(119) 


The  i/’s  in  equation  (119)  are  coefficients  that  define  the  inflow  in  the  same  way  that  the  A’s  define  the 
inflow  in  equation  (118).  However,  the  u's  define  the  inflow  in  the  D  system,  which  is  aligned  with  the 
velocity  in  the  disk  plane.  The  ;v  coefficients  and  the  A  coefficients  are  related  by 


I'D  =  Ao 

•  An 

^Ic  =  ^  +  Ai^C^J 

where  the  side-slip  angle  and  rate  are 


(120) 

(121) 

(122) 

(123) 

(124) 

(125) 


(126) 

(127) 


The  thrust,  roll,  and  pitch  coefficients  in  equation  (1 19)  are  given,  in  terms  of  the  main-rotor  aerodynamic 
forces  and  moments,  by 


(128) 

_  {Sm,,  •  va)  Si3^  +  {Sm^  •  £a)  CfS^ 

(129) 

__  {Sma  ■  va)  Cft  -  (%,  ■  x,,i)  5^, 

(130) 
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The  matrices  for  the  Pitt/Peters  model  are 


Li  = 


1 

Jut 


0 


0 


0 


I  ^  ~  So, 


IStt  /i- 


A//  = 


mQ(ma.ss)  0  0 

0  ^  0 
0  0  5^. 


where 

^  +  •^(A  +  Aq) 

UT  =  \^M^  +  A2 


(131) 


(132) 


(133) 

(134) 

(135) 


The  element  mo  takes  on  two  possible  values.  One  value  corresponds  to  a  zero-inflow  boundary 
condition  at  the  hub  center,  mo(0)  =  jyl;  the  other  value,  mo(l)  =  does  not  enforce  this  condition. 


Tail  Surfaces 

Up  to  five  tail  surfaces  are  available  (the  subscript  representing  each  surface  will  be  omitted  in 
the  following  discussion);  they  are  oriented  as  in  figure  13.  Each  surface  can  be  tipped  away  from  the 
horizontal  by  an  angle  ax,  and  the  pitch  incidence  can  be  adjusted  through  ay  to  define  the  zero  point 
of  the  input  ixs-  The  velocity  of  the  air  relative  to  the  surface  is 

Vs==Vo  +  dTS~^i3  (136) 

where  dxs  is  the  downwash  from  the  main  rotor.  Each  tail  surface  has  its  own  downwash  table  that  is 
a  function  of  the  angles  at  which  the  air  strikes  the  hub  center,  the  magnitude_^of  the  air  velocity,  and 
the  thrust  coefficient  (i.e.,  dxs  =  CxfDWrs  (<^u,^v,Ut,)Ss)-  The  moment  at  the  surface  is  zero, 
but  the  force  at  the  tail  surface  is 

Qs  —l^PA  II  11'^ "■  (<^i) 

+  [—C'l  (Oi)  C'a,  —  C'p  (Or)  5o  J  *3^Qlosa^R  (137) 
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where,  for  forward  flow,  ax  and  S/j  are  defined  by 


ax  =  as 

Sr  =  1 

and,  for  reverse  flow,  they  are  defined  by 


|for-f 


<  <  •! 


=  a  — TT  1  r  TT  ^  ^ 

S«  =  0.8 

=  a  +  7r\.  .  - 

=  0.8  I  fof 2  - 


(138) 


az 

Sr 


TT 

"I 


where 


.  -1  f  -V’i  • 

Qi  =  tan  y— 

V  —  VJj  •  a-5 


(139) 


The  lift  and  drag  coefficients  are  functions  of  oj  and  are  defined  by  specifying  the  break  points  shown  in 
figures  14  and  15.  The  dynamic  pressure  loss  due  to  fuselage  blockage  is  approximated  by  the  functions 


= 1  - 


^  =  -i 


/af  ^  /I3f  1 

\  °Ts  )  \  i'rs  /  J  J 


(140) 

(141) 


When  the  angles  oy  and  are  at  and  6ors»  respectively,  the  loss  is  a  maximum. 

Tall  Rotor 

The  tail-rotor  orientation  is  shown  in  figure  16;  the  Euler  angles  px  and  py  have  the  same  sequence 
as  the  tail-surface  orientation.  The  relative  velocity  of  the  air  is  given  by 


V,.  =  Vq  +  djR  -  ri 


ir 


(142) 


and  the  downwash  cI-^r,  like  the  tail  surface  downwash,  is  interpolated  from  a  table.  From  the  figure, 
the  advance  ratio  and  total  inflow  are 


= 


(Vr  •  x-r)  +  (Vr  •  yj) 


^TR^TR 


^TR  =  I^TR  tan  Or  + 


TR 


—  ^fSTR  4" 


C't- 


TR 


’^'Jy^TR  +  ^TR  ^  '^'J^TR  +  ^TR 

The  wind  orientation  angles  with  respect  to  the  hub  are  defined  as 

-1  I -Vr  •  yT\ 


(143) 

(144) 


/3,.  =  tan’ 


-VV  ■  fr, 


(145) 
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Qr  =  tan“* 


[Vr  •  Xj)  +  (Vr  •  J/rj  J 

The  force  and  moment  for  the  hub  from  the  tail  rotor  are 


ZT 


n2 


(146) 


Qr  =  (-^TR^pr  ~  ^TrSp^  XT  +  {-HtrS/}^  +  YtrCh^)  Vt  -  TjRfblkh  (147) 

Jr  =  [LtrC'ij^  -  MtrS0^)  XT  +  (LtrS0^  +  MtrC0^  vt  +  Qtr^  (148) 

All  components  of  the  force  and  moment  expressions  result  from  the  quasi-static  flapping  model  of 
reference  6,  and  arc  functions  of  the  total  inflow  and  the  blade  harmonic  components.  The  inflow 
comes  from  equation  (144)  and  the  blade  harmonics  result  from 

[/C]a  -  /  =  0  (149) 

where  is  a  vector  containing  the  collective,  the  first  sine,  and  the  first  cosine  components  of  the  blade 
flapping  motion.  The  precise  definitions  of  the  matrices  and  the  rotor  forces  and  moments  are  given  in 
appendix  C  of  reference  6.  In  the  present  report  the  total  inflow  is  defined  such  that  positive  inflow 
is  down  through  the  rotor,  whereas  in  reference  6  total  inflow  is  defined  such  that  positive  inflow  is 
up  through  the  rotor.  Both  equations  (144)  and  (149)  are  algebraic  nonlinear  equations  and  are  solved 
using  the  Newton  method.  Once  acceptable  convergence  is  reached,  the  values  of  ^  and  \tr  can  be 
used  in  the  force  and  moment  equations.  The  loss  factor  in  equation  (147)  is  included  in  order  to 
model  the  interference  effect  between  the  rotor  and  the  surface  to  which  it  is  mounted.  Loss  factor  f^ik 
is  defined  as  follows: 


h\k  "  (T’JZ/oss 


/  ^TR 
xTRbrk 


2 

+  1 


hlk  -  1 


PTR  ^  TRbrk 
UTR  >  TRbrk 


(150) 

(151) 


The  variable  TRi^gg  is  the  percentage  of  thrust  available  at  zero  advance  ratio,  and  TR^rk  is  the  advance 
ratio  at  which  f^ik  becomes  unity  and  remains  unity  at  advance  ratios  above  TR^rk- 


Main<Rotor  Downwash 


The  downwash  of  the  main  rotor  on  the  tail  surfaces  and  tail  rotor  is  calculated  using  the  flat-wake 
model  described  in  detail  in  chapter  2  of  reference  5.  The  wake  model  assumes  a  rotor  in  edgewise 
flow,  as  shown  in  figure  17;  the  coordinate  frame  is  that  of  reference  5.  The  downwash  velocity  is 
determined  using  the  Biot-Savart  law,  assuming  cycloidal  vortex  filaments  shedding  from  the  disk  plane. 
This  flat-wake  assumption  is  considered  valid  in  forward  flight  when  V/RCl  >  1.63v'^.  The  x,  y, 
and  z  components  of  the  velocity  are  calculated  from 


5 

dp 


dp 


x^y, 


(152) 
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where  Atij  results  from 


Avx  =  At)x(xi,yi,Si) 

At>y  =  Auy,  (ii,yi,ii) 

-  + A/(oo,yi,; 

1  fF  r- 

^^2  =  ~2  [A.'(®i,yi,2i)  +  A'(oo,yi,fi 


si)]  +  ('^liyi.si)  +  iV(oo,?yi,ii)]| 
)]  +  [A(xi,yi,xi)  +  i((x,yi,5i)]| 


(153) 

(154) 

(155) 


The  functions  in  equations  (153)  through  (155)  are  pretabulated  and  depend  on  the  nondimensional 
distance  from  the  hub  center  p  and  the  components  of  the  point  of  interest,  through 


il  = 


(156) 


A  typical  example  of  each  of  these  functions  is  shown  in  figure  18  for  yi  =  0.6.  The  functions  need  to 
be  tabulated  only  at  the  positive  values  of  xi,  yi,  and  zi  because  of  the  symmetric  and  antisymmetric 
properties  of  the  functions.  A  summary  of  the  odd  and  even  behavior  of  the  functions  is  shown  in 
table  2  along  with  the  equation  numbers  and  figure  numbers  that  define  the  functions  in  reference  5. 
The  downwash  components  are  expressed  in  terms  of  the  function  T,  which  is  the  average  circulation 
over  a  rotor  revolution  at  a  radial  distance  from  the  center  of  the  rotor.  As  in  appendix  C  of  reference  16, 
this  function  is  approximated  as 


r(?) =ro(i  -  p)T>‘ 

To  = 


5  H  -  4  (ri  -  rf ) 


(157) 

(158) 


where  the  constant  Tq  is  a  function  of  the  thrust  coefficient.  This  relationship  between  the  circulation 
and  the  thrust  coefficient  comes  from  evaluating 


T=  r[;l 

Jr\  [27r  yo  dp 


dp 


where 


/o  dp 

which  results  from  the  Joukowsky  lifting  law. 


(159) 

(160) 


The  wake  model  can  be  used  to  generate  a  downwash  table  for  each  tail  surface  and  tail  rotor.  The 
wake  is  oriented  along  the  relative  air  velocity  vector  in  the  xy  ~  VV  plane,  shown  in  figure  19.  The 
S-frame  is  the  tail-surface  frame,  but  the  tail-rotor  frame  T  can  be  substituted  for  it.  The  velocity  at  the 
hub  is 


Vv  =  Vq  +  -  Uh 


(161) 
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and  the  orientation  angles  are 


a-u  —  tan 


-1 


Vv  •  % 


xh)  +(K;  ■?///)■ 


0v  =  tan 


(162) 


(163) 


For  the  given  tail  surface  (tail  rotor)  fixed  in  the  S  (T)  frame,  downwash  values  /dWts  ttt 

different  combinations  of  av,/3v,  and  ||t^ul|  are  calculated  and  stored  in  a  table.  These  values  are  the 
downwash  velocities  in  the  %  (z^)  direction  divided  by  Cy.  This  division  of  the  downwash  velocities 
by  Ct  is  based  on  equations  (152)  through  (158),  from  w'hich  the  proportionality  of  Ct  to  the  downwash 
velocity  can  be  deduced.  For  the  tail  rotor,  the  downwash  vector  is 

and  for  each  tail  surface  the  downwash  vector  is 

dTS  =  %  (165) 

The  functions  are  linearly  interpolated  from  values  tabulated  from  a  given  set  of  ay,  f3v,  and 

The  points  that  make  up  the  interpolation  tables  are  selected  in  the  following  manner.  The  ||V^||’s 
are  chosen  directly  by  selecting  the  number  and  the  particular  values.  The  ay’s  are  chosen  by 
selecting  a  set  of  points  apr  and  defining  oy  =  app  +  Oyo:  similarly,  the  /3y’s  are  chosen  by  selecting 
a  set  of  points  0pt  and  defining  The  quantities  Oyo  and  /3yo  are  the  orientation  angles 

that  center  the  wake  on  the  given  surface  (i.e.,  when  %  is  parallel  to  f/ij ).  The  points  app  are  chosen 
by  setting  the  boundary  angles  oilObnd  <  ~^DW  <  0  <  ^DW  <  ^UP  <  ^UPbnd' 

points  /3pj  are  chosen  by  setting  the  boundary  angles  PlObnd  <  0  <  0UP  <  (^UPbnd'  These 

boundary  angles  are  defined  for  the  condition  in  which  the  flat  wake  is  centered  on  the  given  surface.  For 
the  a  boundaries,  the  number  of  divisions  in  the  end  regions,  [oilobnd'^^o]  I^UP^^UPbnd^' 
is  selected  through  na^nd’  number  of  divisions  in  the  middle  regions,  [apo, -A/;)Vn 

[^DVVi^UP\>  selected  through  The  value  ^dw  value  and  precludes  evaluation 

of  the  downwash  on  the  wrtke  layer,  where  the  vortex  singularities  are  present.  For  the  0  boundaries, 
the  number  of  divisions  in  the  end  regions,  [Plobnd'I^^o]  »nd  [0[rp, 0uPbnd^'  selected  through 
^/3end’  number  of  divLsions  in  the  middle  regions,  [/?£,o,0]  and  [0,/?t/p],  is  selected  through 

^0mid ' 


Mass  Matrix 


The  equations  of  motion  are  foitned  numerically,  making  the  model  described  in  this  report  implicit. 
Because  of  this,  the  system  mass  matrix  is  not  readily  available  as  it  would  be  if  the  equations  were 
derived  explicitly.  Allowing  for  different  hinge  sequences  precludes  the  explicit  calculation  of  the  mass 
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matrix,  since  each  of  the  many  hinge  sequences  would  require  a  separate  calculation.  To  avoid  this,  the 
mass  matrix  is  implicitly  calculated  with  the  condition  that  no  acceleration-dependent  applied  forces, 
such  as  the  apparent  mass  terms  in  blade-section  aerodynamics,  are  present  in  the  system.  Such  forces 
contribute  to  the  mass  matrix,  making  its  extraction  more  difficult.  It  is  well  known  from  analytical 
dynamics  that  the  form  of  the  kinetic  energy  is 

'^kin  -T2  +  T1  +  To  =  0]  £  +  C{q,  tfq  +  r{q,  t)  (166) 

where  is  the  mass  matrix  of  the  system  and  g  represents  the  system  generalized  coordinates.  In 
the  model  pre.sented, 

rp 

2  ~  P\  1  /^2i '  *  •  1 Cl )  Clf  •  •  •  t  Cni) I  0 1 1 02 1  •  •  •  i  i  i  i  ,  Rnix  i  '  Rmx  i  '0  (167) 

and  the  quadratic  contribution  T2  to  the  kinetic  energy  is  given  by 

Tl  =  +  Tl^otor  +  E  i^'2.bladt)k 

I 

For  the  fuselage,  the  velocity  and  angular  rate  are  given  by 

Wp 

-  ^FU  ^  “  ’^F/I  ^  ®  {  ^  } 
i-lP  ~  '^FiR  +  =  'f'mc<^FII  •  '^FI  |  ^  | 

Rp 

which  is  a  re-expression  of  equation  (22)  using  matrix  notation  instead  of  vector  notation.  The  elements 
of  the  column  matrix  R  are  the  time-derivatives  of  the  three  translational  degrees  of  freedom  describing 
the  motion  of  point  m  on  the  fuselage;  the  elements  of  ^  are  the  time-derivatives  of  the  three  Euler 
angles  of  the  fuselage.  For  the  coordinate  system  F  defining  the  fuselage  motion, 

1  0  ■ 

^F/I=  0  C'fl,  (171) 

_  0  -Sj,  ^ 

The  velocity  and  angular  rate  give  the  T2  contribution  to  the  kinetic  energy  from  the  fuselage 

jP 

Tlfaselagc  =  ^  |  +  W^IfWf]  |  ^  |  (172) 


(169) 

(170) 
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Likewise,  for  the  rotor  hub  mass,  the  matrix  representation  of  the  velocity  and  angular  rate  is 


ai/J;  =  +  W;?///  =  '1'rF'^F/I^  +  =  [^RF^Ffl  •  ^  =  ^R/F 


\Vn 


^  \ 

R 

^  \ 

\  0  J 

=  Tir/jR  +  r[;^i^F/lt  =  \r\2^F/I  i  T./j  i  o’ 


where 


^R/F  = 


Rr 

■  0 
0 
1 


< 

i  '0  J 

(173) 


(174) 


■0  =  ^R/F'4' 


(175) 


Equations  (173)  and  (174)  arise  from  equations  (14),  (18),  and  (63).  As  with  the  fuselage,  the  T2 
contribution  to  the  kinetic  energy  from  the  rotor  hub  is 


'^^roior  2  ] 


f 

T 

<  R 

■0  J 

MRB.jiRji  -H  W^rIrWr 


SL 

R 

0  J 


(176) 


The  same  approach  is  used  to  find  the  kinetic  energy  for  the  fcth  blade.  For  the  Aith  blade  the  angular 
rate  is 

=  ‘J^Rji  +  <J^ii^R  +  iiifl/i  =  Tbf^F/I^  +  T^BR^R/Fi'  +  '^BL^L/R^L/R  +  ^BjL^BIL 


Tbf^F/I  •  0 :  Trr^r/f  :  Tbl^lir  :  u^b/^ 


0 

R 

0 

H/r 

[  ^b/l  j 


(177) 


and  from  equation  (74)  the  velocity  of  the  center  of  mass  becomes 

.(/?)  .(/?)  ,  .(ft) 

Ue  =  Ah  +4e  S^ft// 

.(L)  .(L)  ,  ,(L) 

A/  =  ^L/I 


Ab^  =  4/ ^  +  AfT^B/I  =  [^'h^F/I  =  :  ^'h^R/F  •  ^h^L/R  =  ''')V^'BIL  I  ^ 


(178) 


ftfli. 


0 

R 

0 

'^BfL 
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where 


Mi=TBLfif+f^fTTBL 

M2  =  MiTLfi  +  Tbrv^^ 
M3  =  MiTrf  +  TBF^mh 


Since  no  particular  hinge  sequence  has  been  stipulated,  the  hi^B,  ^^id  4s/ l  generic,  representing 
the  generalized  coordinates  in  the  transformation  indicated  in  the  subscripts.  For  the  Apt,  fplt,  and  pflt 
hinge  sequences,  the  terms  are  deAned  as 


i^L/R  = 


-C, 


\p 


0 


h  =  ^L/RS^L/R 


(179) 


^B/L  = 


^B/L&B/R 


0  1 

L  ^9k-+<j>k  ® 

and  for  the  Ifpt,  Ipft,  and  plft  hinge  sequences,  the  terms  are  deAned  as 


{^t}  = 


(180) 


^L/R  = 


0 

0 

1 


4  =  ^hrh/r 


(181) 


^B/L 


^h^<l>kCip  0 

L  ^h+<t>kC^p  0  J 


\  ^  ,r. 


I  =  ^LIR^B/R 


(182) 


Using  the  velocity  and  angular  rate,  the  T2  contribution  to  the  kinetic  energy  from  the  A:th  blade  becomes 

T 


(^^hlade)  k  O 


Q. 

4 

t/l 


^LjR 

[  ^B/L  J 


[UbRI.Rb^  +  WbJbWbi, 


Q 

4 

V 

^L/R 

[  ^b/l  j 


(183) 


For  the  inAow  model,  the  mass  matrix  associated  with  the  Arst-order  states  xj'  =  [Ao,Ai5,Aic|  can  be 
easily  extracted  using  equation  (119)  and  equations  (123)  through  (125).  The  mass  matrix  for  the  inAow 
model  is 

1  0  0 


=  qLjMi 


I  0  -%  J 


(184) 
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If  there  are  no  acceleration-dependent  applied  forces,  the  total  mass  matrix  of  the  system  can  then  be 
related  to  the  equations  of  motion  by 


Mqq 

0 

0 

Mdi  . 

{ih- 

leq 

ieq 

(£,i,£,^,£j  ^ 

The  matrix  Mqq  is  formed  using  the  T2  terms  in  equation  (168),  by  adding  the  contributions  of  the 
corresponding  fuselage,  rotor  hub,  and  blade  elements.  The  equations  included  in  r^q  are  the  rotor-blade 
equations,  equations  (55)  through  (57)  or  equations  (58)  through  (60);  the  equations  included  in  are 
the  fuselage  equations,  equations  (8)  through  (13)  and  equation  (54).  Finally,  the  equations  included  in 
3^q  are  the  inflow  equations,  equations  (119).  The  reader  will  note  that  these  equations  have  been  cast 
such  that  they  are  identical  to  equations  that  would  result  from  a  Lagrangian  formulation.  Thus,  the 
mass  matrix  calculation  is  compatible  with  the  equations  of  motion  and  allows  the  relationships 


-^4?  - 


-Mdi  = 


dq 

_a 


Leq 


L 


9  J 


=  < 


Uq 

Leq 

[  jC, 


Multiblade  Coordinate  Transformation 


(186) 

(187) 

(188) 


For  rotor  systems  with  three  or  more  blades,  a  multiblade  coordinate  transformation  (MBCT)  is 
available  to  transform  the  rotor  rotating  coordinates  to  nonrotating  coordinates  (ref.  7).  For  a  set  of 
rotating  variables,  each  of  which  represents  the  same  physical  quantity  for  each  blade, 

=  [xi,a;2,a-3, .  ..]  (189) 

The  nonrotating  representation  for  an  even  number  of  blades  is 

x^  =  [xq,  Xj^,  X  jg,  X  ic,  X2s , . . .]  (190) 

The  generic  variable  xq  is  the  collective,  while  xj  is  the  differential,  which  is  not  present  if  the  number 
of  blades  is  odd.  The  Xna  terms  are  the  nth  sine  variables,  while  the  x„c  terms  are  the  nlh  cosine 
variables.  An  example  is  the  flap  degree  of  freedom,  for  which  ii/j;  is  a  vector  whose  components  are 
flapping  angles.  The  transformation  between  variables  is  linear  and  time-varying,  so  the  transformation 
from  rotating  to  nonrotating  systems  is  given  by 


(191) 

(192) 

(193) 


3LN  =  1'nR^R 

i-N  =  Tt\R£.R  +  T^R£.R 

£.N  =  TmrUr  +  '^-T^rLr  +  Tf^R£.R 
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and  the  transformation  from  nonrotating  to  rotating  systems  is 


=  TrnUN 

(194) 

i-R  =  '^RNS.N  +  TrN^N 

(195) 

^R  =  '^RNi.N  +  '^'i'NRi-N  +  '^NRS.N 

(196) 

ITie  form  of  the  transformation  matrix  for  models  with  an  even  number  of  blades  is 


■  1 

-1 

^201 

C'201 

1 

1 

5202 

'^RN  = 

1 

-1 

5^3  C^3 

*  . 

1 

1 

and  its  inverse  is  given  by 

1 

1 

1  1 

-i 

1 

-1  1 

25^, 

25^2 

2503 

'^RN  =  '^NR  = 

1 

2^0, 

2C^ 

2C03 

Tlfj 

232^,2 

-C'201 

(197) 


(198) 


The  V’jfc  terms  are  the  azimuth  positions  of  the  blades: 

^TT 

V'fc  =  +  (fc  -  1) —  A:=  l,2,...,ni,  (199) 

rih 

These  positions  all  have  the  same  time  derivatives.  For  models  with  odd  numbers  of  blades,  the  second 
column  of  equation  (197)  and  the  second  row  of  equation  (198)  are  not  present. 

OPERATIONS 
Basic  Setup 

A  driving  routine  performs  four  basic  operations  on  the  helicopter  mathematical  model:  initial¬ 
ization,  trim,  linearization,  and  integration  (fig.  20).  Initialization  sets  the  basic  data  that  define  the 
helicopter  configuration  and  conveys  to  the  driving  routine  all  the  data  necessary  to  operate  on  the 
model.  The  trim  operation  sets  the  helicopter  configuration  into  some  prescribed  steady  condition.  A 
harmonic  balance  technique  is  used,  which  simultaneously  balances  the  fuselage  forces  and  moments 
and  finds  the  rotor-blade  equilibrium  solution.  Once  the  model  is  trimmed,  two  possibilities  are  avail¬ 
able.  One  is  linearization,  in  which  a  two-point  difference  formula  is  used  to  linearize  the  model  about 
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the  trim  condition.  The  other  possibility  is  time  integration,  which  gives  the  system  time  response, 
initialized  in  trim  condition,  for  given  control  inputs. 

Each  operation  uses  a  subset  of  the  model,  whose  general  structure  is  shown  in  figure  21  and 
represented  in  equation  form  as  follows: 


Leq 

£(X)  =  {  Leg  \  =  -[M{rJ,^,rJ,9,iL) 

3Leq 


r 

f 

X 


^  m  " 
Y=i  V 
X 


(200) 


(201) 

(202) 


Variables  and  controls  are  represented  in  a  single  vector  X.  which  acts  as  the  i^ut  to  the  model,  whereas 
the  equations  and  measured  quantities  are  represented  in  a  single  vector  =  [FJY^],  which  acts 
as  the  output.  Both  the  inputs  and  the  outputs  can  be  categorized  into  three  types:  rotor-blade,  fixed, 
and  augmented.  The  rotor-blade  inputs/outputs  are  items  associated  with  the  rotating  coordinate  system, 
such  as  the  flap  degrees  of  freedom  or  the  blade  pitch;  the  fixed  inputs/outputs  are  items  associated  with 
the  nonrotating  system,  such  as  fuselage  degrees  of  freedom  or  fuselage  acceleration  measurements.  The 
augmented  inputs/outputs  are  associated  with  any  additional  first-order  equations  necessary  to  complete 
the  model.  The  structure  in  figure  21  allows  the  use  of  rotating  or  nonrotating  inputs  and  outputs. 
The  variables  are  transformed  by  the  multiblade  coordinate  transformation;  the  choice  of  the  coordinate 
system  depends  on  the  particular  operation.  This  basic  setup  is  general  and  the  operations  can  work  on 
any  model  set  up  in  this  manner. 

For  the  model  described  in  the  Helicopter  Equations  of  Motion  section,  the  rotor-blade  degrees  of 
freedom  in  the  rotating  system  are 

iF  =  (203) 

When  using  nonrotating  coordinates,  the  vector  t  is  defined  according  to  equation  (190)  for  each  degree 
of  freedom.  For  a  model  with  an  even  number  of  blades  the  rotor  degrees  of  freedom  in  the  nonrotating 
system  are 

i-  “  /^Ifli/^lc)  •  •  •  iCOiCd)  ClsiClc)  ■  •  • )  00 1  0c/  )01s)01c)’  •  ■]  (204) 


For  a  model  with  an  odd  number  of  blades,  the  differential  coordinates  are  not  present.  Using  nonrotating 
coordinates  for  the  input  variables  requires  that  the  flag  ISNROT  =  1,  which  transforms  the  nonrotating 
coordinates  to  the  rotating  system  via  equations  (194)  through  (196). 


The  rotor-blade  controls  in  the  rotating  system  are 


(205) 
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and  as  with  the  rotor-blade  variables,  the  function  of  the  flag  IINROT  is  to  transform  9.  to  the  rotating 
system  if  these  inputs  are  given  in  the  nonrotating  system.  For  a  model  with  an  even  number  of  blades, 
the  rotor-blade  controls  in  the  nonrotating  system  are 


c)  •  •  • )  ^0)  ^ls>  ^Ic)  ‘  '  •  )  ^Ot  ^Isi  ^1 


Cl 
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For  a  model  with  an  odd  number  of  blades,  the  differential  coordinates  are  not  present.  The  fixed 
variables  are  the  fuselage  degrees  of  freedom  and  the  rotor  azimuth  angle 

^  dz,  Rrtix  1  Rmyi  R-irii  i  '*/’]  (207) 

and  the  augmented  variables  are  the  three  dynamic  inflow  coefficients 

£  “  ['^Ot  <^lc]  (208) 


The  fixed  controls  are 
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The  center  block  of  figure  21  is  the  model,  described  in  the  previous  section,  that  calculates  the  helicopter 
equations  of  motion.  The  equations  in  teq  are  the  rotor-blade  equations,  equations  (55)  through  (57) 
or  equations  (58)  through  (60),  and  the  equations  in  are  the  fuselage  equations,  equations  (8) 
through  (13)  and  equation  (54).  The  equations  in  are  the  inflow  equations,  equations  (119).  The 
vector  £  is  defined  by  equations  (40)  through  (47),  and  jn  is  defined  by  equation  (105)  and  equations  (107) 
through  (112).  The  vector  w.  is  in  the  rotating  frame  but  can  be  transformed  to  the  nonrotating  frame 
by  using  the  flag  lONROT. 


Initialization 

Initialization  performs  two  basic  functions,  one  of  which  is  to  set  up  any  data  not  given  in  X 
that  define  the  model.  These  data  are  the  physical  parameters  defining  the  model  that  was  previously 
described;  they  are  summarized  in  table  3.  The  main-rotor  downwash  tables  are  either  generated  and 
stored  in  files  or  read  in  from  previously  generated  files.  The  second  function  of  initialization  is  to 
transfer  basic  system  data  to  the  driver  so  that  subsequent  operations  can  be  carried  out  on  the  model. 
These  data,  summarized  in  table  4,  are  primarily  size  data  of  the  full  model.  Included  in  the  table  are 
the  current  values  of  these  data  for  the  model  described  in  the  Helicopter  Equations  of  Motion  section 
of  this  report. 


Trim 

The  trim  operation  places  the  helicopter  in  the  desired  flight  confi‘?uration  by  using  a  harmonic 
balance  technique  that  is  described  in  reference  8,  chapter  3.  The  method  casts  the  helicopter  equations  in 
an  algebraic  form  by  defining  the  rotor-blade  variables  as  a  Fourier  series.  The  coefficients  of  this  series 
expansion,  along  with  the  fixed  and  augmented  state  variables  and  inputs,  are  adjusted  to  force  the  rotor- 
blade,  fuselage,  and  augmented  state  equations  to  zero,  The  rotor-blade  equations  are  also  expressed  as 
a  Fourier  series,  so  that  forcing  the  rotor  blade  equations  to  zero  means  forcing  the  coefficients  of  this 
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series  expansion  to  zero.  It  is  assumed  that  the  helicopter  main  rotor  has  a  constant  rotation  rate  ijj  =  Cl 
and  that  the  rotor-blade  controls  are  expressed  in  the  nonrotating  frame  (i.e.,  IINROT  =  1).  All  blades 
are  assumed  identical,  and  information  from  only  one  blade,  expressed  in  terms  of  rotating  coordinates 
(i.e.,  ISNROT  =  0),  is  needed  in  the  trim.  An  IMSL  nonlinear  algebraic  equation  solver,  ZSPOW,  is 
used  to  find  the  solution  to  the  following  equation  (ref.  17); 

£(Z)  =  Q  (210) 


The  vector  of  inputs,  ]/,  is  defined  as 
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which  is  composed  of  fixed,  augmented,  and  rotor-blade  inputs.  (Matrices  are  stacked  by  columns 
when  represented  in  column  vectors,)  The  blade  variables  are  represented  as  coefficients  of  the  Fourier 
expansion. 


r[k  +  nij{j  -  1)]  =  rcoihj)  + 


^^TRM 
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rco(2nJ)  cos(nV'A;)  +  rco(2n  +  l,i)  sin(nipk) 


(212) 


The  variable  k  is  the  blade  number,  and  j  is  the  blade  degree  of  freedom  (e.g.,  j  =  2  is  the  lead-lag; 
see  cq.  (203)).  For  the  model  of  this  report,  the  coefficients  are  represented  in  matrix  form  as 
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The  elements  of  V,  are  adjusted  by  the  equation  solver  after  the  user  chooses  determining  indices.  For 
example,  the  choice  of  nj^  =  2  with  if{\)  —  2  and  if{2)  =  3  directs  the  solver  to  adjust  angles  By  and 
6z  to  find  the  trim.  Given  the  information  in  Y_,  the  input  vector  X  can  be  formed  and  the  model  can 
be  evaluated.  The  resulting  equations  can  then  be  used  to  form  £,  which  is  defined  as 
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These  equation  are  the  Fourier  coefficients  of  the  first  blade  equation,  where 

^^trm 

reqik  +  TibU -\)]=rfo{lJ)+  Y1  ^/o(2n,i)cos(ni/)fc) +r/o(2?i+ l,7)sin(7ti/)fc)  (216) 

n=l 

The  choice  of  determines  the  number  of  rotor  coefficients  that  are  adjusted  to  force  the  same 

number  of  coefficients  of  the  blade  equations  to  zero.  The  fixed  and  augmented  equations  are  averaged 
over  one  rotor  revolution  in  equation  (214)  to  zero  only  the  constant  harmonics  of  these  equations.  A 
summary  of  the  inputs  that  control  the  trim  procedure  is  shown  in  table  5.  Indices  and  their  number  (in 
parentheses)  are  selected  to  choose  the  adjustable  variables  and  the  equations  that  are  to  be  satisfied.  If 
an  adjustable  variable  is  set  to  a  value,  the  value  acts  as  an  initial  condition  to  the  solver.  If  the  variable 
is  not  adjusted  by  the  solver,  it  remains  at  the  set  value  throughout  the  trim.  An  example  of  the  use  of 
this  feature  would  be  to  trim  the  model  at  a  constant  climb  rate. 


Integration 

The  tine  response  of  the  system  shown  in  figure  22  to  given  inputs  Up  can  be  obtained  by  numerical 
integration.  Since  it  is  difficult  to  generalize  helicopter  control/actuator  models,  these  systems  have 
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been  broadly  defined  as  first-order  equations  that  can  be  defined  by  the  user.  The  P  block  in  figure  22 
represents  precompensation,  and  the  A  block  is  any  actuator/mixer  dynamics;  feedback  dynamics  are 
represented  in  the  C  block  and  sensor  dynamics  are  represented  in  the  S  block,  which  has  the  restriction 
that  the  output  of  S  cannot  explicitly  depend  on  Ug.  The  helicopter  model  is  represented  in  block  G  with 
nonrotating  rotor  inputs  (IINROT  =  1)  and  rotating-blade  degrees  of  freedom  and  outputs  (ISNROT  =  0, 
lONROT  =  0).  By  default,  P  and  A  are  identity  feedthrough  and  C  and  S  are  zero,  giving  the  response 
of  the  open-loop  helicopter  model  G.  Any  subset  of  the  model  represented  by  equation  (200)  can  be 
integrated,  where  the  state  of  G  is  defined  by 
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As  in  the  the  trim  procedure,  the  selection  of  the  variables  and  equations  is  done  by  choosing  indices. 
If  the  rotor  degree  of  freedom  is  not  selected,  a  constant  rotor  rate  is  used,  by  default,  in  the  integration 
of  the  model  (i.e.,  'ip  =  fl).  The  inputs  are  defined  as 


S.{K{n) +nb[j0{m)  -  \]}  n  =  1, 2, . . ,  ,nt;  m  =  1,2, , , 
\  n[iu(m)]  m  =  l,2,..,,n4 

The  feedback  outputs  are 

'  ui{n  +  ni[jw{m)  -  1]}  n  =  l,2,...,n;,;  m  =  1,2,.,. 

VG  =  S  m=  l,2,...,nfj^ 

,X[ix(^)]  m=l,2,.,.,nx 

and  the  auxilary  outputs  are 

{  Ul{n  +  niy\jii,{m)  -  1]}  n  =  1,2, .. .  in  =  1.2, , ,, 
|/[t^(m)]  m=  1,2,...,%- 

X  ['jx(^)]  m  =  1,2, , , .  ,nj^ 

The  equations  integrated  are 
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From  equation  (200)  the  form  of  the  subsystem  that  is  integrated  is 
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where 
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Two  methods  are  available  for  extracting  the  mass  matrix  Mq.  One  method  uses  equation  (222) 
directly,  in  which  it  is  clear  that  the  jth  column  of  Mq  is  Eq,  with  da  equal  to  zero  except  for  the 
jth  element  set  to  unity,  minus  If  is  of  dimension  =  ni,nj^  +  rnj.  -I-  nj,, 

one  evaluation  of  the  model  is  necessary  to  get  Zg(Q,  aic,iiG!),  and  n^Q  additional  evaluations  of  the 
model  are  necessary  to  get  Mq.  This  approach  always  works;  however,  it  requires  many  evaluations 
of  the  model  and  thus  is  inefficient.  The  other  method  is  to  use  the  corresponding  elements  from  the 
mass  matrix,  which  is  calculated  in  and  passed  from  the  model,  to  form  the  submatrix  Mq-  Once  Mq 
is  available,  the  equation 
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can  be  found;  this  is  the  expected  form  of  block  G  in  figure  22. 
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If  the  helicopter  model  is  trimmed,  the  inputs  and  outputs  of  G  are  available  at  the  initial  time.  The 
initial  state  xq^,  initial  input  and  initial  output  j/Gq  can  be  used  to  extract  the  initial  conditions 
of  the  blocks  A,  C,  P,  and  S  (fig.  22)  as  follows.  It  is  assumed  that  =  ojg  =  a;p  =  ^5  =  0  and 
it  is  required  that  nuQ  =  nup>  which  gives  an  equal  number  of  equations  and  unknowns.  With  these 
conditions,  a  nonlinear  equation  (j;)  =  0  can  be  set  up,  where  ^  and  n  are  defined  as  in  figure  23. 
The  variables  on  the  left-hand  side  of  the  figure  are  adjusted  until  the  variables  on  the  right-hand  side 
are  forced  to  zero.  The  values  at  which  this  occurs  are  the  initial  conditions  of  the  system  of  figure  22. 
The  solution  is  extracted  with  the  modified  Newton  method  used  in  the  trim  procedure. 


Wich  the  initial  conditions  available,  the  entire  system  of  figure  22  is  integrated  according  to  the  flow 
diagram  in  figure  24.  The  S  block  is  evaluated  to  get  its  output  yg  using  Ug,  which  is  not  necessarily 
the  the  correct  value  at  the  given  time.  However,  the  output  of  S  is  specified  to  be  independent  of 
H5,  so  the  output  value  is  correct.  With  this  input  to  block  C,  the  values  on  the  right-hand  side  of  the 
figure  can  be  evaluated  from  the  given  information  on  the  left-hand  side  of  the  figure.  With  this  setup, 
and  the  input  up  defined  by  the  user,  the  integration  is  carried  out  using  the  predictor-corrector  method 
in  the  subroutine  DVERK  from  the  IMSL  package  (lef.  17).  A  summary  of  the  integration  inputs  is 
shown  in  table  6.  As  in  the  trim,  the  variables  to  be  integrated  are  selected  through  indices.  The  input 
is  Up  =  Upo  +  ^Upy  where  UjpQ  is  determined  in  the  trim.  The  term  AZp  is  the  input  perturbation 
about  this  trim. 


Linearization 

After  trim,  it  is  also  possible  to  linearize  a  subset  of  the  helicopter  model  (c.g.,  block  G  of  fig.  22) 
in  the  nonrotating  system  (ISNROT  =  1,  IINROT  =  1,  and  lONROT  =  1).  Linearization  is  done  about 
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the  trim  solution  using  the  two-point  finite-difference  formula, 


Z  (^i  +  ~  L  (®j  ~ 

dxj  ~  2Axj 


(225) 


The  partial-derivative  matrix  of  the  helicopter  model  is  defined  in  figure  25;  the  rows  and  columns  are 
selected  in  the  same  way  as  in  the  integration  process.  For  any  given  time,  the  partial-derivative  matrix 
can  be  calculated  and  the  linear  system  can  be  evaluated  according  to  the  following  equation: 
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These  integrals  are  evaluated  numerically  using  a  Gaussian  quadrature  procedure.  The  periodic  model 
matrix  is  then  expressible  as 
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Since  nonrotating  coordinates  are  used,  the  helicopter  model  matrix  will  be  constant  when  the  model  is 
in  hover.  In  forward  flight,  the  model  matrix  has  periodic  coefficients,  which  become  more  important  as 
the  forward  flight  speed  is  increased.  The  fundamental  frequency  w j  depends  on  the  number  of  blades 
in  the  helicopter  system.  If  nj  is  odd  then  Uf  =  ni,Cl,  whereas  if  it  is  even,  uif  =  ^12.  A  summary  of 
the  linearization  procedure  inputs  is  given  in  table  7. 


RESULTS 


The  model  results  are  compared  with  other  experimental  and  analytical  results  as  a  preliminary 
validation.  The  first  comparison  is  with  blade-equilibrium  responses  in  trim,  of  a  flap-lag-torsion  rotor 
model  used  in  reference  8.  Unlike  the  model  described  in  the  present  report,  the  equations  in  reference  8 
were  symbolically  generated  out  to  their  final  fonn.  An  omission  in  the  reverse-flow  region  was 
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discovered  in  the  model  of  reference  8:  the  model  did  not  account  for  the  sign  change  in  in  its 
aerodynamic  formulation  (see  fig.  11).  This  error  is  introduced  into  the  model  described  herein  so  the 
comparison  would  be  valid.  The  parameters  of  the  model,  presented  in  table  8,  are  of  a  four-bladed, 
soft  in-plane  rotor  with  a  lead-lag  frequency  of  0.7/rev.  The  rotor  trim  is  calculated  with  six  harmonics 
{nfiTRM  ~  degree  of  freedom,  and  the  fuselage  mass  is  chosen  to  give  a  weight 

coefficient  of  0.005.  The  flap,  lag,  and  torsion  equations  are  satisfied  for  the  blade  equilibrium  while 
simultaneously  satisfying  the  roll-moment,  pitch-moment,  ;c-force,  and  ::-force  equations.  The  trimmed 
blade-equilibrium  response  in  flap,  lag,  and  torsion  is  shown  in  figure  26.  Three  response  curves  are 
presented  for  each  blade  degree  of  freedom:  (1)  a  curve  from  the  model  from  reference  8,  which  has 
an  erroneous  reverse-flow  formulation;  (2)  a  curve  from  the  model  of  the  present  report,  which  uses  the 
same  reverse-flow  formulation  as  that  used  in  reference  8;  and  (3)  a  curve  from  the  model  of  the  present 
report,  which  uses  what  is  believed  to  be  the  correct  interpretation  of  the  aerodynamic  formulation 
in  the  reverse-flow  region.  The  two  models  with  the  incorrect  reverse-flow  interpretation  show  close 
agreement,  which  serves  to  verify  the  analysis.  The  alternative  reverse-flow  interpretation  greatly  affects 
the  torsion  response  of  the  system,  although  its  effect  on  the  flap  and  lag  responses  is  still  very  small. 

The  results  from  the  linearized  helicopter  model  are  compared  with  experimental  hover  results  from 
reference  9.  The  experimental  setup  consisted  of  a  three-bladed  rotor  mounted  on  a  mast  that  allowed 
pitch  and  roll  motions.  Each  blade  had  metal  plates  with  notches,  which  acted  as  concentrated  hinge 
springs  to  closely  simulate  the  flap-lag  rigid-blade  approximation.  Torsion  springs  on  the  body  mount  of 
the  experimental  rotor  were  adjusted  to  simulate  an  air-resonance  condition.  The  parameters  necessary 
to  model  this  system  are  given  in  table  9.  Measurements  were  made  of  the  damping  and  frequencies 
of  the  system  modes,  and  it  is  to  these  data  that  a  coaelation  is  made.  Figure  27  shows  the  modal 
frequencies  and  the  lead-lag  regressing  damping  plotted  against  the  rotor  rotation  rate,  with  zero  blade 
pitch.  The  simulated  air-resonance  condition  occurs  when  the  body  roll-mode  frequency  nears  that  of 
the  lead-lag  regressing  mode,  which  is  near  a  rotor  rotational  rate  of  700  rpm.  Excellent  correlation  with 
the  experimental  data  is  seen  in  this  condition.  The  body  pitch  and  roll  damping  at  zero  blade  pitch  is 
shown  in  figure  28.  Good  correlation  between  experiment  and  theory  is  also  seen  here,  although  the 
model  tends  to  overpredict  the  damping  at  the  higher  rotor  rates.  The  lead-lag  regressing  damping  at 
9°  blade  pitch  is  shown  in  figure  29;  again,  the  model  does  a  good  job  of  predicting  the  damping  in  the 
resonance  condition.  Detailed  cross-plots  of  the  lead-lag  regressing  damping  with  respect  to  the  blade 
pitch  are  shown  in  figure  30.  Good  agreement  is  seen  between  experiment  and  theory,  although  the 
damping  tends  to  full  short  at  high  and  negative  blade-pitch  angles  at  the  higher  rotor  rates.  Figure  31 
shows  cross-plots  of  the  body  damping  at  650  rpm;  good  agreement  in  pitch  is  seen  at  positive  blade- 
pitch  angles.  The  body  roll  damping  tends  to  be  overpredicted,  getting  worse  at  higher  biade-pitch 
angles. 

The  main-rotor  wake  model  used  to  calculate  the  downwash  on  the  empennage  surfaces  is  tested 
by  correlating  it  with  experimental  results  from  reference  5,  page  77.  A  correlation  is  done  in  the 
reference  and  is  also  done  in  the  present  report  by  using  the  flat-wake  model  described  in  the  Main 
Rotor  Downwash  subsection.  Unlike  the  model  used  in  reference  5,  the  model  described  herein  assumes 
a  circulation  distribution.  In  reference  5,  a  blade-element  approach  is  used  to  find  the  circulation 
distribution.  The  integral  in  equation  (152)  is  evaluated  by  using  a  Gaussian  quadrature  integration 
method  along  with  linearly  interpolated  values  from  the  pretabulated  functions,  Figure  32  shows  the 
vertical  nondimensional  downwash  at  various  azimuths  at  0,17?  below  the  rotor.  The  experimental  rotor. 
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with  10°  blade  twist,  is  immersed  in  edgewise  flow  at  Ct  —  0.006.  Using  the  flat-wake  model  gives 
good  correlation  results;  the  asymmetry  of  the  downwash  at  the  90°  and  270°  azimuth  positions  is 
captured. 

With  the  rotoi  and  downwash  models  validated,  it  is  now  possible  to  validate  the  full  helicopter 
model.  The  model  is  configured  to  simulate  a  UH-60A  helicopter  in  order  to  compare  model  responses 
with  the  flight-test  data  that  are  available  in  references  10  and  18.  Reference  10  presents  a  validation 
of  the  real-time  model  presented  in  reference  11,  which  is  also  known  as  the  GEN  HEL  model.  Ref¬ 
erence  18  presents  a  validation  of  an  updated  version  of  the  GEN  HEL  model,  which  is  referred  to  as 
the  Ames  GEN  HEL  model.  In  the  model  of  the  UH-60A  used  in  this  report,  all  automatic  controls 
are  off  (SAS,  PBA,  etc.).  The  fuselage  aerodynamic  representation  in  the  model  was  generated  from 
curve-fitted  equations  of  the  static  aerodynamic  characteristics.  The  exact  equations  are  given  in  ref¬ 
erence  19,  pages  2-3,  and  are  used  in  place  of  the  equivalent-drag-area  fuselage  aerodynamic  model. 
The  lag  damper  is  modeled  according  to  reference  11,  page  5.1.10,  and  this  representation  is  substituted 
for  the  spring/damper  model.  The  control  mixer  is  modeled  according  to  reference  11,  pages  5.5.14  to 
5.5.17,  without  the  servo  models.  The  tail-rotor  parameters  are  extracted  from  reference  20,  where  a 
detailed  description  of  the  UH-60A  tail  rotor  is  given.  References  11  and  12  are  used  to  extract  other 
physical  parameters  of  the  helicopter.  Table  10  summarizes  of  all  of  the  inputs  to  the  model.  The  data 
followed  by  a  question  mark  are  specified  for  each  of  the  following  correlations.  From  table  10  it  can 
be  seen  that  the  model  has  a  flap-and-lag  main  rotor  and  a  six-degree-of-freedom  fuselage  at  constant 
rpm.  The  dynamic  inflow  model  is  included,  and  the  collective  pitch,  sine  pitch,  cosine  pitch,  and 
tail-rotor  collective  pitch  act  as  the  inputs  to  the  helicopter  model.  No  feedback  outputs  are  selected, 
and  accelerometer  outputs  are  taken  from  the  fuselage  at  the  positions  corresponding  to  those  on  the 
flight-test  vehicle. 

Trim  results  calculated  with  and  without  the  effect  of  the  main-rotor  wake  in  straight  and  level 
flight  are  compared  with  experimental  results  in  figure  33.  The  trim  data  were  calculated  at  25-knot 
increments  with  the  horizontal  tail  surface  at  the  following  settings:  40°  for  0,  25,  and  50  knots;  11.3° 
for  75  knots;  and  4°  for  100,  125,  and  150  knots.  The  vehicle  was  trimmed  by  satisfying  all  six  fuselage 
force  and  moment  equations,  the  flap  and  lag  rotor  equations,  and  the  inflow  equations  (see  table  10). 
Between  0  and  50  knots,  the  angle  Oz  was  fixed,  while  the  angle  0^  was  adjusted  (i.e.,  if  =  1,2  in 
table  10).  Between  50  and  150  knots,  the  angle  6x  was  fixed,  while  the  angle  Oz  was  adjusted  (i.e., 
•i/  =  2,3  in  table  10).  The  collective,  longitudinal,  lateral,  and  pedal  inputs  are  presented  in  percent 
of  total  motion;  the  total  motion  is  10  in.  for  the  stick  inputs  and  5.38  in,  for  the  pedals.  The  model 
tends  to  deviate  significantly  from  the  experimental  data  above  120  knots,  where  aerodynamic  effects 
such  as  dynamic  stall  and  radial  flow,  which  are  not  accounted  for  in  this  model,  can  become  important. 
Below  120  knots,  good  correlation  is  seen  and  the  results  are  similar  to  those  presented  in  reference  18, 
p.  28.  The  lateral  and  longitudinal  stick  positions  are  predicted  more  accurately  by  the  Ames  GEN  HEL 
model,  but  the  pedal  position  is  closer  for  the  model  described  herein.  Overall,  the  results  from  the  Ames 
GEN  HEL  model  are  slightly  better,  probably  because  that  model  used  tabular  data  for  blade-section 
aerodynamics,  as  well  as  other  refinements  particular  to  the  UH-60A  configuration.  The  flat-wake  model 
enhances  correlation  in  all  but  the  collective  and  lateral  stick  positions.  The  collective  stick  position  is 
hardly  changed,  whereas  the  lateral  stick  position  is  degraded  by  1/4  in.  at  150  knots. 

Flight  test  time-histories  from  reference  10  are  compared  with  the  responses  of  the  UH-60A  model, 
at  1  knot  (figs.  34-37)  and  at  100  knots  (figs.  38-41).  Four  sets  of  time-respon.ses  are  presented  at  each 
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flight  speed;  they  correspond  to  1  in.  adjustments  to  the  collective,  longitudinal,  lateral,  and  pedal 
inputs.  The  inputs  to  the  model  are  those  from  the  flight  data,  referenced  about  the  trim  of  the  model. 
The  trim  values  generated  by  the  model  are  indicated  by  points  c,  e,  a,  and  p,  corresponding  to  the 
collective,  longitudinal,  lateral,  and  pedal  inputs,  respectively.  The  angular  rates  of  the  fuselage  (roll, 
pitch,  and  yaw)  and  the  accelerometer  outputs  from  the  fuselage  (longitudinal,  lateral,  and  vertical) 
are  also  presented  for  each  case.  The  accelerometer  position  corresponds  to  the  approximate  center  of 
gravity  of  the  helicopter;  the  precise  center-of-gravity  location  is  given  as  POSFUS  in  table  10.  Each 
plot  has  three  curves:  the  test  data,  the  model  response,  and  the  linearized  model  response.  The  test 
data  are  from  reference  10  and  the  test  numbers  are  identical  to  the  numbers  used  in  the  reference.  The 
model  responses  come  from  integration  of  the  full  nonlinear  system,  and  the  linearized  model  responses 
come  from  integration  of  the  linearized  model.  The  responses  of  the  linearized  model  are  added  to  the 
trim  values  to  give  the  total  responses  presented.  Since  the  helicopter  is  in  straight  and  level  flight, 
the  angular  rates  are  zero  at  trim,  and  the  linearized  responses  are  identical  to  the  total  responses.  The 
accelerometer  responses  contain  gravity  effects,  and  the  nonzero  trim  values  of  the  accelerometers  must 
be  added  to  the  linearized  accelerometer  outputs. 

At  1  knot  (figs.  34-37),  the  on-axis  angular  rates  agree  well  with  the  test  results.  The  off-axis 
responses  of  the  model  do  not  correlate  as  well,  with  the  pitch  axis  being  the  most  troublesome  by 
diverging  from  the  test  data  in  tests  212,  203,  and  209.  However,  this  is  no  worse  than  the  results 
from  the  GEN  HEL  model,  which  is  refined  specifically  to  the  UH-60A  configuration  and  includes  an 
engine  model  with  a  rotor  rpm  degree  of  freedom  and  actuator  dynamics.  At  100  knots  (figs.  38-41)  the 
on-axis  responses  have  fair  agreement  with  experimental  results,  and  again,  these  results  are  consistent 
with  results  from  the  GEN  HEL  model.  Other  tests  were  run  at  60  knots  and  140  knots,  and  the  results 
show  a  strong  similarity  to  the  results  from  the  GEN  HEL  model.  Another  source  for  comparison  is 
reference  16,  in  which  the  results  of  a  linearized  model  are  compared  with  the  same  flight  data.  The 
results  from  reference  16  are  also  similar  to  results  obtained  with  the  model  of  this  report.  Taking  all 
these  results  together,  it  is  concluded  that  the  present  model  gives  results  that  are  consistent  with  those 
obtained  from  the  GEN  HEL  program  and  the  linearized  model  of  reference  16. 

A  final  comparison  is  made  of  the  linearized  model  with  the  eight-state  quasi-static  model  generated 
from  the  Ames  GEN  HEL  program.  The  linearized  GEN  HEL  model  is  generated  with  Q  =  21  rad/sec, 
ST,n  =  STc  =  STa  =  29.583  ft,  WL,n  =  WLc  =  WLa  =  20.683  ft,  Mp  =  492.13  slugs,  iaero  =  0, 
and  cry  ■—  [0, 0, 0]  rad.  The  linearized  model  of  this  report  is  generated  and  arranged  according  to 


y  =  [^r  M  (232) 

The  vector  x.r  contains  all  rotor  states  and  inflow  states,  and  the  vector  Xf  contains  all  fuselage  states. 
The  inputs  are  the  longitudinal,  collective,  lateral,  and  pedal  positions  {uj  =  [(i,.,  be,  ba,  (^i/)]^).  The  out¬ 
puts  ^  are  cho.sen  to  correspond  to  the  states  of  the  GEN  HEL  model  {{^  =  [u,  tr,  q,  A9,  v,  p,  A0,  r,  A0]). 
For  comparison  with  the  GEN  HEL  model,  the  derivatives  of  the  rotor  and  inflow  states  iv  are  set  to 
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zero,  yielding 
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(238) 

The  state  variables  used  to  describe  the  fuselage  motions  are  not  precisely  those  used  in  the  GEN  HEL 
model;  however,  the  information  needed  to  make  the  transformation  to  the  GEN  HEL  states  is  found  in 
the  equation 
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The  quantities  u,  t;,  and  w  are  the  velocities  along  the  fuselage,  and  p,  and  r  are  the  body  angular 
rates.  The  quantities  A0,  A0,  and  Aijj  arc  the  pitch,  roll,  and  yaw  angular  perturbations.  With  this 
choice  of  output  quantities,  it  is  clear  that  Bred  =  0.  and  thus  equation  (239)  holds.  Equation  (239)  is 
then  used  to  transform  the  model  to  the  final  form 
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where  has  been  removed  from  the  system. 


(240) 


This  operation  was  carried  out  on  a  linearized  model  at  1  knot;  the  resulting  F  and  G  matrices 
are  presented  in  table  11  along  with  the  system  matrices  from  the  Ames  GEN  HEL  program.  The 
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eigenvalues  and  normalized  eigenvectors  are  snown  in  figure  42  for  both  models.  The  calculations  are 
based  on  the  F  matrix  represented  in  dimensions  of  feet,  seconds  and  degrees  so  as  to  scale  the  angular 
rates  to  an  order  of  magnitude  comparable  to  the  velocities.  Real  eigenvalues  have  real  eigenvectors 
(figs.  42(d)  and  42(e));  complex  eigenvalues  have  complex  eigenvectors  composed  of  a  real  part  and 
an  imaginary  part.  The  eigenvector  corresponding  to  the  complex  conjugate  of  the  eigenvalue  is  the 
complex  conjugate  of  the  eigenvector.  The  greatest  difference  between  models  arises  in  the  higher 
frequency  modes,  shown  in  figures  42(a)  and  42(b).  The  model  of  this  report  tends  to  introduce  more 
coupling  between  roll  and  pitch  than  the  Ames  GEN  HEL  model.  The  mode  in  figure  42(c)  is  dominated 
by  w  and  r  in  both  models,  although  the  content  differs.  In  figures  42(d)  and  42(e)  the  two  modes  are 
very  close,  although  the  damping  is  much  higher  in  the  model  of  this  report.  Differences  in  the  damping 
and  frequencies  of  the  eigenvalues  are  also  noted;  they  are  likely  caused  by  the  different  inflow  models 
used  in  the  two  analyses.  Selected  Bode  plots  of  the  system  are  shown  in  figure  43.  The  6c/ w  plot 
shows  good  agreement  and  the  6p/r  gain  agrees  well  with  the  Ames  GEN  HEL  model.  The  360°  phase 
shift  in  the  6p/ r  plot  is  due  to  a  zero  in  the  right  half  of  the  s-plane  that  is  not  unstable  in  the  modci  of 
this  report.  The  sharp  peak  in  the  6a /p  plot  is  due  to  the  light  damping  that  is  predicted  by  the  Ames 
GEN  HEL  model. 

The  entire  process  is  repeated  at  100  knots,  resulting  in  the  F  and  G  matrices  given  in  table  12.  The 
eigenvalues  and  eigenvectors  are  presented  in  figure  44;  the  eigenvectors  agree  more  closely  between 
models  in  this  case  than  in  the  1-knot  case.  The  on-axis  Bode  plots  of  the  100-knot  case  are  given  in 
figure  45.  The  two  models  do  not  agree  closely,  but  the  general  trends  in  the  models  are  similar. 

CONCLUDING  REMARKS 

A  single-rotor  helicopter  mathematical  model  for  the  design  of  high-bandwidth  flight  control  sys¬ 
tems  has  been  described.  The  flap,  lag,  and  torsion  rotor  model  allows  the  approximate  modeling  of 
hingeless  rotor  systems,  and  the  various  hinge  sequences  that  are  available  allow  the  modeling  of  many 
different  articulated  rotors.  Quasi-steady  Greenberg  aerodynamics  are  used  in  conjunction  with  the 
nonlinear  Pitt/Peters  dynamic  inflow  model  to  calculate  the  blade  aerodynamic  forces,  An  rprn  degree 
ol  freedom  of  the  main  rotor  is  available,  along  with  a  six-degree-of-freedom  rigid  fuselage.  The  tail 
surfaces  have  main-rotor  downwash  effects  that  arise  from  the  flat-wake  simulation.  Operations  on  the 
mode!  include  trim,  linearization,  and  time-integration,  which  have  been  formulated  so  that  they  may 
be  applied  to  any  consistently  cast  rotorcraft  mathematical  model.  This  features  allows  for  modification 
of  the  helicopter  model  vdthout  requiring  significant  changes  in  the  operations.  Preliminaiy  validation 
of  the  model  showed  reasonable  correlation  with  experimental  and  analytical  results. 

Further  validation  of  the  model  is  necessary.  The  higher  order  linearized  models  need  to  be 
compared  to  experimental  results,  and  flap-lag-torsion  experimental  data  are  needed  to  validate  the 
torsional  effects  of  the  rotor  model.  Also,  additional  work  is  necessary  to  extend  the  range  of  valid  flight 
conditions,  which  are  currently  limited  to  low-speed,  lightly  loaded  helicopter  systems.  The  limiting 
factor  is  the  blade-section  aerodynamic  model,  which  should  be  improved  to  account  for  dynamic 
stall,  radial-flow  effects,  and  compressibility  effects.  Finally,  the  linearization  process  is  limited  to  the 
helicopter  model  itself.  A  useful  additional  operation  would  be  the  linearization  of  the  remaining  blocks 
of  the  system  and  the  combination  of  these  linearized  blocks  into  a  single  set  of  system  matrices. 
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Table  1 .  Rotation  sequences  of  blade  coordinate  systems 
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Table  2.  Function  odd/even  properties 
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Table  :b  Summary  of  helicopter  properties 
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Variable  Value 
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Valu 

Rotor 
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Table  9,  Flap-lag  rotor/fuselage  data 


Vanable 


Rotor/fuselage 
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U.209  kg 
0.811  m 
?  rad/sec 


0  0.010069 
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0,0851  m 
9.81  m/sec” 
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0.726  m 
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10 
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1,2761  kg/m^ 

0.02618  rad 
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0.0079 
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6,691  kg  m^/rad  sec^ 

30.659  kg  m^/rad  sec^ 

0.003538  kg  m“/rad  sec 
0.007574  kg  m^/rad  sec 
68.03  kg  m^/rad  sec^ 

104.3  kg  m^/rad  sec- 
0,081 17  kg  m^/rad  sec 
0.4200  kg  nr^/rad  sec 
'0.183  0  Ol 

0  0.633  0  kg  m- 

0  0  0  j 

20.83  kg 
0.241  m 


Variable 
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kg  m- 


^(1.3) 


MPTLIN 


Value 


Trim 


2 

30 

2 

0,01  m/sec 


do  rad 
1 
2 
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[1.2] 

[1.2.3] 

[3] 

[2.3] 

[1.2] 

Linearization 

0 

10 

2 

[1.2] 

2 

[1.2] 

2 


[1.2] 

3 

[1.2,3] 

[10“'.  10“^]  rad/sec 


[10  10  ^1  rad/sec 

[10-^  10“^]  rad 
[10“^  10“^.  10“^] 


Table  10.  UH-60A  configuration  data 
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0  0  0  ■ 
0  573.42  0 

0  0  573.42 
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-0.05236  rad 
?  ft 
?  ft 
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20 

26.83  ft 

[20,60,100,140,180,220]  ft/'sec 
3 


[1,2,3,4.5,6,7] 


Linearization 
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0-1  DM  0.1066  0.4<J50  -0.8769 
0.0595  -0.1700  1.4011  -0.4950 


Table  12.  System  matrices  from  GEN  HEL  program  and  model  of  this  report  at  100  knots 
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(^AMES  <JEN  HEL 


FORWARD 


^  ’’im  “  ’‘l 
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^Fl]  =  ^SrSg^Cg^  -  Cg^Sg^  Sg^Sg^Sg,  +  Cg^Cg^  Sg^Cg^^ 

.  ^h^dyCe,  +  Sg^Sg^  Cg^Sg^jSg^  -  Sg^Cg,  Cg^X'Qy  _ 

V//  =  “  ^zSgy)if  +  {^yCg^  +  B:Sg^Cg^  )ijf  +  {9:Cg^Cg^^  -  BySg^  )i 
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Figure  1.  Fuselage-fixed  coordinate  frame  (F). 
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Rgure  2.  Free-body  diagram  of  fuselage 
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Figure  8.  Lag-flap-pitch-torsion  hinge  sequence. 
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Figure  11.  Blade-section  aerodynamics. 
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Figure  12.  Actuator-disk  coordinate  frame  (A), 
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Figure  13,  Tail-surface  orientation  and  aerodynamic  forces, 
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Figure  20,  Overview  of  operations. 
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Figure  26.  Blade  responses  from  the  present  model  compared  with  those  from  a  flap-lag-torsion  rotor 
in  trim. 
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Figure  27.  Modal  frequencies  and  lead-lag  regressing  damping  data  from  the  present  model  compared 
with  those  from  a  flap-lag  rotor/fuselage  system  at  zero  blade  pitch, 
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Figure  28.  Body  pitch  and  roll  damping  data  from  the  present  model  compared  with  those  from  a 
flap-lag  rotor/fuselage  system  at  zero  blade  pitch. 
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Figure  29.  Lead-lag  regressing  damping  data  from  the  present  model  compared  with  those  from  a 
flap-lag  rotor/fuselage  system  at  9*  blade  pitch. 
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Figure  30.  Lead-lag  regressing  damping  data  from  the  present  model  compared  with  those  from  a 
flap-lag  rotor/fuselage  system  at  various  rotor  speeds  and  blade  pitches, 
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Figure  31.  Body  pitch  and  roll  damping  data  from  the  present  model  compared  with  those  from  a 
flap-lag  rotor/fuselage  system  at  650  rpm. 
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NON-DIMENSIONAL  RADIUS 

Figure  32.  Downwash  from  the  wake  model  compared  with  experimental  results. 
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Figure  40,  Concluded. 
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Figure  41.  UH-60A  model  responses  with  1 -in. -right  pedal  at  100  knots. 
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Figure  42.  Eigenvalues  and  eigenvectors  of  this  report’s  model  and  the  Arres  GEN  HEL  model  at 
1  knot. 
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Figure  43*  On-axis  Bode  pic 
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Figure  45.  On-axis  Bode  plots  of  present  model  and  Ames  GEN  HEL  model  at  100  knots. 
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